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ABSTRACT 


The wind driven surface of a large anechoic tank was used 
to study forward scattered underwater sound. A 10 kHz saw- 
tooth signal was used to drive an omnidirectional source to 
insonify the rough surface. Direct and reflected path sound 
Signals to a point hydrophone were separated by appropriate 
gating and the reflected path signal was analyzed for up to 
2 minutes at 10 and 20 msec. intervals, uSing digital FFT 
spectral analysis. Graphs of sound pressure level vs. time 
for every 10 kHz from 30 kHz to 120 kHz (Surface roughnesses 
between 1.1 and 10.5) are shown. Correlation of the scattered 
sound pressure with frequency was analyzed. Evidence is pre- 
sented to support the hypothesis that, for large roughnesses, 
the instantaneous scattered sound amplitudes are harmonic 
functions of sound frequency; that is, there is a consistent 
frequency separation between instantaneous amplitude maxima 


at one frequency and minima at another frequency. 
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Pee RO wUe Ll LON 


For a number of years, acousticians have studied the 
fluctuations of underwater acoustic signals when scattered 
from a statistically rough interface, typically either the 
surface or bottom. This effect has been described as "build- 
ing and fading" of the signal; an analogy has frequently been 
drawn between this effect and the twinkling of stars. The 
mechanism of this scattering has been the object of many 
theories; References 1 through 5 refer to the basic develop- 
ments in this area. The theoretical developments have common 
shortcomings: in order to solve the integral equation gov- 
erning the scattered sound field, it is necessary to make 
many assumptions involving the boundary conditions, the value 
of the normal derivative on the surface, the spatial statis- 
tics of the insonified area; most analyses deal with high or 
low frequency limits; shadowing and multiple scattering 1s 
usually ignored in all but the most complex and sophisticated 
treatments. These comments concerning the shortcomings of 
strictly theoretical approaches to the scattering problem have 
been made in order that the need for supplemental experimental 
research be noted. A mixture of theoretical insight and con- 
trolled experimental procedure is sometimes more fruitful 
than purely theoretical studies with arbitrary assumptions. 


Such is the philosophy adopted for this research. 
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Tae OBJECTIVES 

The objective of this research is to investigate the 
effects that variation of frequency has on the sound scattered 
from a statistically rough, fluctuating, wind driven water 
surface, and to determine if a predictable relationship exists 
between the intensities of various harmonics of a given sound 
source when scattered from such a surface. The basic con- 
cepts of this experimental effort are: 


1. Insonification of the fluctuating surface with many 
frequencies simultaneously. 


2. Isolation of surface reflected sound from other paths 
through the use of gated signals. 
BIA INTERMEDIATE GOALS 
The following intermediate goals are presented as logical 
steps in attainment of the objective described above: 


1. Perform preliminary studies of the wind driven surface, 
and monofrequency sound scattered from such a surface. 


2. Develop the electronic circuitry to transmit a wide 
range of frequencies simultaneously, receive and process 
these Signals. 


3. Develop and test an acoustic source capable of trans- 
mitting over a wide range of frequencies, with the stip- 
ulation that frequencies of interest be far enough above 
the 'noise' to permit meaningful signal processing. 


4. Analyze the processed data to see if the stated objec- 
tive can be met. 


Inar 





iM eotoe AKC FACTLITIES 


A. OCEAN ACOUSTIC WAVE FACILITY (OAWF) 

The Ocean Acoustic Wave Facility is a combination of a 
wind generated wave tunnel and an anechoic acoustic tank. 

The wind tunnel is approximately 17 meters long, 1.2 meters 
wide and 1.2 meters deep. By varying the number of fans and/ 
or the distance between the water surface and the top of the 
tunnel, a wide range of ocean surfaces can be obtained. 
Figures 1 and 2 (from Tourville) show the essential features 
of the OAWF. 

The wave generation tunnel empties into a water filled 
tank with dimensions of 3 meters on each side. The wind gen- 
erated waves are dissipated on a constructed beach to prevent 
reflected waves from interfering with wind generated waves. 
The tank is constructed of 4x4 redwood posts, placed with 
corners facing outward. The use of redwood provides sound 
absorption, as shown in Reference 8; however, scattering of 
sound from these surfaces did present some experimental diffi- 
culties, as pointed out in a later section. Equipment was 
positioned within the tanks by use of movable racks as shown 
in Figure 3. The OAWF models the real ocean on a scale on the 


eeoom ©: 150 to 1 up to sea state 3. 


5B. DATA ACQUISITION AND PROCESSING CAPABILITIES 
Data acquisition and processing were accomplished by 
utilizing a digital computer system composed of three primary 


components each of which is interfaced to provide high speed 
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analog to digital conversion, digital processing, and data 
printout and recording. The design was developed by the 
Special Projects Section of the Naval Air Development Center 
iMmecmyuNCcLION with Pinkerton Computer Consultants, Inc., of 
Warminster, Pennsylvania. The three components are pictured 
in Figure 4 and described below. 
meeentcerdata Model 70 Computer 
This minicomputer is a digital design that is FORTRAN 
programmable with a 64 thousand byte memory. In addition to 
actual core memory, data that have been stored on digital 
cassettes can be read into the computer for processing. 
2. Enoenix Analog to Digital Converters, Model ADC 712 
Two ADC 712 analog to digital converters may be used 
separately of Simultaneously. Each converter is a high speed 
device capable of encoding £ 10 volt input signals in digital 
formeanad providing an accuracy of 0.005 volts in 20 volts. 
The maximum sampling rate is 320,000 samples per second; this 
rate was routinely used. 
3. Texas Instruments Silent Electronic Data Terminal 
Model 733 
The TI 733 consists of a keyboard used as a program- 
Ming input/output control device, a printer, and a playback/ 
record section used in conjunction with digital cassettes. 
The overall system facilitates rapid, accurate processing of 
any desired type of analog electrical signal and was Used 
mainly for statistical and frequency domain analysis uSing 
Standard FFT algorithms. Further descriptions are contained 


in References 6 and 7. 


ILS 
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Figure 4. View of OAWF 
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4. IBM-360 Interface 
In addition to the digital system just described, 
another processing capability consisted of the TI 733 termi- 
nal and taped digital data interfaced via phone link to the 
IBM-360/67 CP/CMS system. This capability was used to pro- 
vide data rearrangement, card punch capability and graphical 


Satput (with the aid of a CAL-COMP DPlOEECT) . 


C. STANDARD EQUIPMENT LIST 
Much of the equipment used was standard off the shelf 
scientific measurement equipment; a list of standard equip- 


ment referred to throughout the text is given below: 


Abbreviation ue besera peLon 

Scope Tectronix Model 545B 4-channel 
oscilloscope 

HP 465A Amp Hewlett-Packard Voltage Amplifier 
Model 465A 

HP 466A Amp Hewlett-Packard Voltage Amplifier 
Model 466A 

HP 467A Amp Hewlett-Packard Power Amplifier 
Model 467A 

HP 218 Hewlett-Packard Model 218 Digital 
Delay Generator 

HP 219 Hewlett-Packard Model 219 Dual 
Pulse Unit 

PAR 113 Amp Princeton Applied Research Amplifier 
Model 113 

hae 3322 Filter Krohn-Hite Model 3322 Freq. Filter 

KH 3342 Filter Krohn-Hite Model 3342 Freq. Filter 

fe 3350 Filter Krohn-Hite Model 3350 Freq. Filter 

Freq. Counter Hewlett-Packard 5223L Electronic 


Frequency Counter 


al 





DMM 


GR-1192B 


GR=-1217C 


GR-1396 


Impedance Bridge 


Wavetek 114 
Wavetek 136 


Wavetek 144 


Fluke Model 8000A Digital Multi Meter 


General Radio Electronic Frequency 
and Period Counter 


General Radio Model 1217C Unit Pulse 
Generator 


General Radio Model 1396A Tone 
Burst Generator 


General Radio Model 1650A Impedance 
Bridge 


Wavetek Model 114 Function Generator 
Wavetek Model 136 Function Generator 


Wavetek Model 144 Function Generator 


Sound sources, receivers, wave probes and other miscellaneous 


equipment will be described in the text, as necessary. 
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fii. BAPERIMENTAL PROCEDURES 


The experimental procedures followed in the course of 
this research can be grouped into four general categories: 
(1) Calibration of the Surface Wave Probe and eee eninseson 
of wave spectra, wave heights and surface roughnesses; 

(2) Continuous wave sound scatter measurements; (3) Sound 
Source testing and selection; (4) Surface Reflected Pulsed 


Sound Wave analysis. 


A. WAVE PROBE CALIBRATION, WAVE SPECTRA AND WAVE HEIGHT 

MEASUREMENTS 

In two previous works [Refs. 6 and 7], wave height fre- 
quency spectra of wind generated water waves were determined. 
Confirmation of these previous studies was considered EOubDe 
an excellent starting point of this current research and 
would afford valuable insight into the use of the specialized 
measuring equipment and into the use and programming of the 
digital minicomputer. It was intended to extend these 
investigations over those previously performed by studying 
waves generated by all combinations of fans and also to 
develop a program to give actual RMS values of wave height. 

1. Wave Probe Calibration 

The wave height probe is a linear capacitative probe 

developed by Stanford University and consists basically of a 
thin wire whose total capacitance is determined by the water 


depth to which the probe is immersed. The probe is connected 
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to the impedance bridge and the remainder of the circuitry 
as shown in Figure 5. The impedance bridge used had a con- 
venient feature which provided a 1 kHz output Signal modulated 
by the value of the input capacitance. As the probe depth is 
varied the final signal (which is pure D.C. for a smooth 
surface) provides a linear output proportional to the probe 
depth. The probe was mounted on a rigid frame, but was able 
to be varied in depth by means of a rotating knob-screw 
arrangement. A graduated depth scale was attached to the frame 
with an accuracy of about 0.02 cm. The probe was then cali- 
brated as shown in Figure 6. Initial attempts at calibration 
were not successful due to lack of repeatability. This was 
corrected by slight modifications to the filtering network 
and the final calibration curve was satisfactory. 
2. Wave Spectra Measurements 

Upon completion of the wave height probe calibration, 
measurement of wave height spectra was undertaken using the 
Circuitry of Figure 5 as shown previously. The time varying 
Signal proportional to instantaneous wave height was input 
into one A/D channel of the digital computer for analysis 
using the SKIP-1C program (all programs used are described 
briefly in Appendix A). The computer A/D Sampling rate was 
set at 20 Hz; 50 runs of 256 (2°) Samples each were taken at 
this sample rate. Previous spectral measurements had indi- 
cated that most of the spectral components were below 4 Hz 
so that the Nyquist criterion was satisfied in all cases. 


Spectral resolution for this particular data gathering 
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arrangement is 0.078 Hz. For each of the 50 data sets, 
Fast Fourier Transform (FFT) analysis was performed and the 
output spectral density for each run was computed, normalized 
and averaged. The output spectral densities of surfaces for 
all fan combinations of one to five fans were computed and 
plotted. The graphs of output power spectral densities vs. 
frequencies are shown in Figure 7. This set of curves shows 
the expected spectral shifts as the wave heights and wave 
periods are increased by the use of more fans. The general 
results were in agreement with Perkins and Tourville, hence 
further data processing simply for the purpose of curve smooth- 
ing was not performed. As a further check on the computer 
accuracy and program, signals were also input to a PAR Fourier 
Analyzer and the spectral output vs. frequency were plotted 
on an XY recorder. Spectra produced by this alternate means 
were consistent with those analyzed by digital means. A typ- 
ical example of this alternate spectral curve is shown in 
pigure 3. 
me oaverage Wave Height Determination 

After the wave height probe calibration curve was 
obtained, the next task to be undertaken was accurate deter- 
mination of average wave heights for all fan combinations. 
This was accomplished by using the circuitry of Figure 5 
once again, in conjunction with a computer program entitled 
CAL Ol. This program sampled the wave height output voltage, 
performed root-mean-square averaging of the total number of 


A/D data points, and compared this with a reference sin wave 
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voltage set exactly to 1.000 volts. (This comparison elimi- 
nated the need for A/D voltage calibration each time the 
program was used.) This program was extensively tested with 
Signals of known RMS values to ensure accuracy prior to 
experimentally determining the wave heights. Once the pro- 
gram was known to be working properly, RMS voltage values 

for all fan combinations were determined, and converted to 
wave heights using the linear calibration curve. These values 


are shown in Table I. Once the average wave heights were 


Table I 


RMS Wave Heights vs. Number of Fans 


No. of Fans AC (RMS) Volts RMS Wave Height (Cm) 
i Oz0475 0.447 
2 Or Gaks 0 0.546 
5 0.0244 0.740 
4 0.0343 1.041 
5 0.0416 e267 


obtained for each fan combination, it was then a simple matter 


to compute surface roughness values using the definition: 


kk 
= a= (29cos 4) 


where ae is defined as the roughness ee Oo 1s the RMS 
wave height, A is the wavelength of the sound source and 4 
is the angle of incidence. The roughness values were com- 
puted using the geometry of Figure 21; values are contained 


mm Table II. 
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Table If 


Roughness Parameters 


Freq (kHz) Roughness Yq 
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B. CONTINUOUS WAVE SOUND SCATTER MEASUREMENTS 


It was decided to initially analyze the scattered sound 
spectra ofa continuous 20 kHz source. The source and re- 
ceiver setup and associated circuitry are shown in Figure 9 
and are Similar to that of Perkins so that comparison of re- 
sults could be made. The source used was a 2" ceramic 
Sphere; the receiver was a Standard LC-10 hydrophone. The 


recelver was positioned in a Lloyds' mirror maximum signal 
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meoure 9.Sound Processing Circuitry 
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location. The demodulated scattered sound envelope was input 
into the digital computer and analyzed with the SKIP-1C pro- 
gram. The normalized scattered sound spectral densities 

for 1 to 5 fans are shown in Figures 10-14. Each spectral 
density curve 1s an average of 50 FFT's, each consisting of 
256 sample points, sampled at a rate of 128 samples/sec. 

Once again results were consistent with previous experimen- 
ters' results. Oscilloscope photos of typical modulated and 


demodulated sound signals are shown in Figures 15 and 16. 


eee SOUND SOURCE SELECTION AND TESTING 

As discussed in the introductory chapter, in order to 
insonify an area simultaneously with a wide range of frequen- 
cies, it was necessary to test various sources to see which 
could provide the maximum transmitting frequency response in 
the range of interest, nominally 5 to 160 kHz. 

meeleansducer Construction and Testing 

A significant number of potential transducers were 
constructed and tested; for the sake of brevity only a brief 
description of each will be given along with the reason for 
rejection. 

a. Two-inch ceramic sphere, hollow interior: this 
transducer was utilized in the initial sound scattering 
Studies at 20 kHz. The transducer did not have a significant 
Sutpuc below 20 kHz and significant ‘ringing’ occurred at 
50 kHz which would have made the sphere unsuitable for short 


separation gated pulses. 
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Figure 10. Surface Modulated Sound Spectrum,1 Fan. 
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Figure ll. Surface Modulated Sound SSpeceuum, 7 Fans 
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Figure 12. Surface Modulated Sound Spectrum, 3 Fans. 
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FIGURE 13. Surface Modulated Sound Spectrum,4 Fans. 
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Figure 14. Surface Modulated Sound Spectrum, 5 Fans. 
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Figure 15. Surface Modulated Sound Signal 





Figure 16. Surface Demodulated Sound Signal. 
(Figure 16 is not the demodulation 
of Figure 15,but was taken at a 
different time and is a typical 
demodulated signal.) 
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b. One-inch ceramic sphere, hollow interior: same 
comments as in a., extensive ringing around 75 kHz. 

c. One-inch ceramic sphere, oil filled: same prob- 
lems as previously discussed. Oil was introduced into the 
cavity to damp out the ringing phenomenon but it only shifted 
the resonance values. 

d. One-inch ceramic sphere, oil filled with steel 
wool damping. No major improvement noted here. 

e. One-inch ceramic sphere, packed with heavy 
powdered metal. Same comments as c. apply. 

f. Four cylindrical mylar transducers were con- 
structed, 1" and 2" in diameter, 10" and 20" in length respec- 
tively. The transducers had an aluminum core, thin mylar 
covering, with the entire assembly coated with neoprene water- 
proofing compound. When tested, all four of these transducers 
exhibited too low a signal level to be useful. 

g. LC-10 ceramic hydrophone. This hydrophone ex- 
hibited satisfactory signal level over the entire frequency 
range; due to its small size and relative omnidirectional 
features, this standard receiving hydrophone was utilized as 
an acceptable source throughout the remainder of the experi- 
ment. Actual characteristics of this transducer will be 


shown and discussed in subsequent sections. 


D. SURFACE REFLECTED PULSED WAVE ANALYSIS 
This section will detail the experimental procedure used 


to analyze the surface reflected sound as a function of 
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frequency and time. Prior to studying the behavior of sound 
in the anechoic tank, a study of the sawtooth oscillator was 
initiated as described below. 
1. Wavetek 144 Sawtooth Waveform Analysis 

As previously mentioned, a Wavetek 144 function gen- 
erator was modified to give an asymmetrical sawtooth wave- 
with a slope ratio of 44:1. This asymmetrical output was 
coupled directly into the A/D input channel of the digital 
computer. The THCDB program was utilized for Fourier analyses 
of this input signal. Each entire analysis of the 10 kHz 
sawtooth waveform consisted of 49 sample blocks, each block 
consisting of 256 sample points taken at a sample rate of 
320 kHz, giving an available frequency resolution of 1250 Hz. 
The sawtooth waveform was analyzed in both a continuous and 
gated mode. Examination of the output data (mean spectral 
density values, and standard deviation vs. 5 kHz frequency 
increments) revealed an unacceptably high value in the higher 
frequency ranges. Values of standard deviation in dB re 
1 ubar for frequencies over 100 kHz were greater than 3 dB 
in all cases. In order to make subsequent experiments mean- 
ingful, a method had to be found to decrease the standard 
deviation to less than 1 dB for all frequencies to be 
examined. 

The first study performed was to compare the Wavetek 
144 spectral analyses with similar results from other signal 
generators. As can be seen in Table III, three different 


models of function generators gave similar results, showing 


the increased instability of harmonics with frequency. 
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Various methods were tried in order to overcome this 
difficulty; amplitude variation was investigated but had no 
Significant effect. The ultimate solution to this problem 
proved to be band pass filtering of the output of the Wavetek 
prior to analysis. An extensive study of the effects of var- 
iation of filter limits was performed; results showed that 
filtering was effective in reducing the standard deviation 
of the output signal to within 0.5 dB or less. Further 
studies of this nature were performed with the direct and 
reflected sound paths in the wave tank and will be subse- 
quently discussed. 

Ba circuit Arrangement for Sound Analysis 

The basic circuit arrangement for all future dis- 
cussions on sound analysis is shown in Figure 17. A brief 
discussion of the basic function of each component is neces- 
Sary to understand the processing of the various sound sig- 
nals. The Wavetek 144 function generator controls the basic 
repetition frequency of the entire transmissions circuit. 

The GR 1217 Pulse Generator controls the actual length of the 
Sawtooth waveform, and the positive going portion of its 
apse initiates the timing sequence of the HP 218/219. The 
HP 218/219 Digital Delay Generator, when operated in the 
‘common!’ mode has as its output two pulses of variable ampli- 
tude and width, each of which can be set to a pre-adjusted 
time delay in microseconds. This component also has a syn- 
chronizing output for use in oscilloscope 'synch' input. The 


time delayed pulses are used as the timing or gating signals 
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for the GR-1394 Tone Burst generator. This component, in 
turn, allows the passage of the 320 kHz sampling frequency 
to the computer only as the gating pulses permit. 

The GR-1217 also controls the sawtooth waveform 
going to the LC-10 source. The receiving circuit consists 
of an LC-10 hydrophone, 80 dB of amplification, and a band- 
pass filter. The output Signal of the filter goes to the 
A/D input channel for computer analysis. The series of 
oscilloscope photographs shown in Figures 18-20 present the 
various waveforms and their interrelation. In examination 
of these figures, the appearance of the waveforms should be 
noted. In Figures 19 and 20, the direct path and the surface 
reflected path can easily be identified; however in both 
cases, additional scattering is present following the surface 
reflected wave, consisting of a mixture of acouStic energy 
scattered from the side and side-surface of the wave tank. 
It was the presence of this additional scattered wave which 
determined the minimum separation (about 7 msec.) between 
successive sound pulses and also determined the maximum length 
the surface reflected pulse could attain (about 0.6 msec.) 


without overlap with the scattering. 
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IV. EXPERIMENTAL RESULTS AND ANALYSIS 


A. STATIC ACOUSTIC WAVEFORM INVESTIGATION 

Upon completion of the testing and spectral analysis of 
the direct output of the Wavetek 144 oscillator, testing of 
the direct path and smooth surface reflected acoustic sig- 
nals was initiated. Although some variation of geometry was 
initially utilized, the geometry of Figure 21 was used for 
all analysis unless otherwise noted. The circuit arrange- 
ment of Figure 17, shown in the previous section, applies. 

A study of the spectral components of a 10 kHz sawtooth 
waveform transmitted via the two different paths of interest 
was conducted. Although filtering of the signal was per- 
formed, the old problem of high standard deviation among the 
components was noted but could not be reduced to acceptable 
levels; many of the standard deviations were on the order of 
2 dB and this was felt to be unacceptable for the supposedly 
Stable direct and reflected path signals. Further study re- 
vealed that both the PAR-113 and KH-3550 had voltage limita- 
tions on their output, which if exceeded, resulted in clip- 
ping of the output. The PAR-113 started to clip at = 6 volts 
and the KH-3550 at an output of = 8 volts. By reducing the 
magnitude of the input signal, the amplified received sound 
could be kept within these bounds; when the signals were kept 
within these limits standard deviations were greatly reduced 
and were acceptable (Table IV). Due to the inherent limita- 


tions of the PAR-113, it was replaced by a HP-465 amplifier 
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Figure 21. Source and Receiver Geometry. 
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which had no limitations near the range of voltage values 
encountered. Particular care had to be taken when analyzing 
mae cluctuating signals of the rough surface to ensure the 
clipping point of the KH-3550 was not approached. During 
Emus phase Of testing, the filter limits of 20-130 kHz were 
‘found to produce the most stable signals with the least 
overall standard deviation of the direct and reflected 
waveforms. 

The spectral output of the direct and surface reflec- 
ted sound paths, as analyzed by the THCDB program are shown 
in Table IV below. These values are representative for the 
smooth surface case. Typical noise values, obtained by dis- 
connecting the source, and analyzing the received signal with 
the THCDB program are also shown. The results of this table 
are shown graphically in Figure 22. Each mean value and 
standard deviation were computed over a 50 block average, with 
each block consisting of 128 sample points sampled at 320 kHz. 
The slightly higher standard deviation for the reflected path 
at 120 kHz is thought to be due to a minor surface rippling 
due to slight air movement over the wave tank. Note that most 
Signal to noise ratios for the reflected wave are on the order 


of 30 dB or higher. 


B. ROUGH SURFACE SOUND SCATTER MEASUREMENT 

Surface sound scatter measurements were now attempted for 
all combinations of 1 to 5 fans. As before, the circuit 
arrangement of Figure 17 was used. Gated sound pulses were 


set to a 20 msec. interval between pulses; the A/D sampling 
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Figure 22. Reflected Waveform Spectrum,Smooth Surface. 
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pulses were positioned in time to be directly coincident with 
reception of the surface scattered sound pulses. Each of the 
50 blocks consisted of 128 A/D samples taken at a rate of 

320 kHz; thus each sample block took 0.4 msec. During this 
small time interval, it was assumed that the rough wave sur- 
face is frozen in space over the sampling period. This assump- 
tion is justified based on examination of the wave frequency 
Spectra previously shown geese the peak energy 1S at a period 
of approximately 300 msec. 

The output sound pressure levels vs. frequency for each of 
the time separated reflected pulses was computed using the 
THCDB program, and stored on a separate digital tape cassette 
for each fan combination. This tape was, in turn, processed 
by another program, SORTF, which sorted the SPL vs. time by 
frequency. At this point, for each 5 kHz increment of frequency 
there existed a set of SPL data over 50 points, each point 
separated by 20 msec. in time, encompassing a total time span 
of 1 sec. The data of the SORTF program was also stored on 
digital tape cassettes. Utilizing a phone link with the IBM 
360 computer complex, the SORTF data was card punched and sub- 
sequently plotted using the IBM 360 plotting routines and the 
CAL COMP plotter. 

ime intetiakl Results 

Typical results from the previously described process 
are shown in Figures 23 and 24. The particular curves shown 
are for relatively low surface roughness values. Study of the 


whole series of curves obtained for the 1-5 fan cases indicated 
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the existence of a possible problem. The pseudo periodicity 
of the SPL curves was on the order of 12-15 Hz; this was 
approximately three times that expected for the lower rough- 
ness values. A complete review of experimental and computa- 
tional procedures was initiated to seek the answer to this 
problem. A reexamination of the THCDB program showed that 
after the A/D data for one sample had been stored, Fourier 
analysis and all remaining computational processes were being 
performed prior to the storage of the next A/D data set. [In 
effect, what resulted was a skipping of the timed computer 
sampling orders until the computational steps were performed, 
which resulted in timed samples being taken only every 60 to 
80 msec. vice the 20 msec. separation which was expected. As 
a result of this finding, it was decided to modify the exist- 


a 
ing programs to eliminate the problem. 


2. Program Modifications 

In order to get data at the precise time desired, a 
number of program modifications were made as described below: 

a. Raw A/D data were taken and stored on a digital 
cassette using an existing program, PERKI1A. In order to ensure 
the data were being taken at the exact time desired, testing 
was performed on 10 kHz triangle wave with A/D data being 
taken with a 5 msec. separation. The resulting A/D data were 
manually plotted and fitted exactly to the predicted triangle 
form. The PERKIA program was subsequently modified to be able 
to take up to 100 data blocks of 128 points each. The modi- 


fied program was called PRKMOD. 


Sel 





b. The THCDB program was modified to be able to 
accept input data from a cassette and to be able to compute 
and store on tape the output spectral analyses of up to 100 
blocks in time. The modified program was tested on a static 
reflected wave from the smooth surface. This was compared 
with the original program and nearly identical results were 
obtained; the modified program was labeled THCDB(1). 

c. The SORTF program was modified to accept up to 
100 data blocks. 

See tinal Experimental Data 

With the program modifications described above final 
experimental data taking was commenced. The first set of data 
taken was for fan combinations of 1 to 5, pulse and sample 
time separation of 20 msec., for a total duration of 1 sec. 
of data. The resulting graphs of this data were plotted with 
the aid of the CAL COMP plotter and a redrawn overlay showing 
curves for selected fan combinations is given in Figures 25, 
25a, 26. The amplitudes shown are relative and placed on each 
overlay such that the mean values of individual curves are a 
fixed distance (2") apart; these mean values are shown on the 
overlays as ‘tic marks' on the ordinate axis. It will be noted 
that the 5 fan combination is missing; during 5 fan data taking, 
a Slow undulation of the direct path sound signal was noted; 
this undulation was not present in the 1 to 3 fan cases but 
was also noted in the 4 fan case after about 10 minutes at this 
fan combination. It 1s surmised that this variation in the 


direct path sound was due either to relative motion of the 
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source/receiver induced by wave slap on the mounting rods or 
by some type of oscillating current within the wave tank itself. 
For these reasons, 5 fan data were no longer taken. 

The examination of the first data set indicated that 
the time resolution in the upper frequencies might not be fine 
enough; consequently an additional set of data with a 10 msec. 
time separation was taken encompassing a total time of 0.75 
seconds. These data were plotted as before and overlays shown 
fame dgures 27, 28, and 29. 

It was also felt that the total time display of data 
might not be long enough to show changes in the SPL which were 
taking place; therefore a last set of experimental data was 
taken, utilizing a 20 msec. pulse Separation. but over nearly 
100 points for a total time history on the order of 2 seconds. 
These data were plotted and overlays shown in Figures 30, 3l, 
amc 32. 

Analysis of the entire data array is contained in sub- 
paragraph 5 of this section. Individual plots showing actual 
sound pressure levels and accurate time scales are available 
in the Ocean Physics Laboratory, Naval Postgraduate School. 

4. Post Experimental Data Processing 

After completion of the graphical presentation of the 
three varied data sets, and after being faced with the some- 
what overwhelming mass of this graphical data, it was felt that 
further data processing might lead to some simplification of 
analysis. The form that this further processing was to take 


was that of cross-correlation; a method was desired by which 
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a SPL vs. time curve at one frequency could be compared with 
that af another frequency and reduced to a single number 
representing the relative similarity of the two curves. Thus, 
instead of correlating in time or space, a correlation was em- 
ployed using the algorithm below; it should be noted that this 
is a rather special type of cross-correlation wherein products 
of pressures at one frequency are compared with pressures at 


another frequency but summed over time. 





T2 
> Ca (Me ur (PK (2) >= uF.) 
cr, 5) i F 2 
Corr(F,F,) = 2 7 7 
2 
: (Pp (4) = uF) OS (Pp, (i) = uF)? 
é Py il Es 2 
1=T) i=T, 
The parameters of this algorithm are defined as follows: 
Py = frequency fora given correlation curve 
= frequency to which Fy is being compared or correlated 
Po, (1) is the sound pressure (in microbars) at the ith time 
1 interval for frequency Fy 
Pa (i) is the sound pressure (in microbars) at the ith time 
2 interval for frequency F 
UP, = mean value of Pa over the total time interval of concern 
1 
UF, = mean value of Pe over tht total time interval of concern 
2: 
T = starting point in time 
T, —wmencingepoint in time. 


Each correlation curve for a given Fy consists of a set of ll 
correlation values (30 kHz to 130 kHz, taken at 10 kHz inter- 
vals). It should be noted that the curves are normalized, 


1.e., when Fy = Poe Corr(F)F.) = l. 
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The actual processing of the data was relatively simple 
Since all data cards had been previously punched. A new com- 
puter program was developed for the IBM-360 to carry out the 
correlation algorithm and provide the appropriate output. A 
copy of this program is contained in Appendix A. 

In addition to the correlating over the entire number 
of time points available, correlation was also performed over 
the halves or over individual fourths of some curves for com- 
parison purposes. Table V describes those correlations which 
were computed and plotted. Complete presentation of all the 
correlation curves produced (ll curves for each case listed 
in Table V) is not considered feasible just from a numerical 
standpoint but representative results are labeled and shown in 


Figures 33 to 39. 
ris 


>. Analysis of Results 


Examination of Figures 25-31 gives insight into the 
qualitative aspects of the effects of frequency change on the 
scattering of sound. The relative periodicity, corresponding 
to the periodicity of the wave surface, is present for the 
eeeively low roughness, smoother surface cases. As the fre- 
quency, and hence the roughness, increases this periodicity 
disappears. It is assumed that the incoherent component of 
the sound is growing more predominant and is responsible for 
this effect. Actual measurements to separate coherent and 
incoherent sound components were not made. For the 1 and 2 
fan overlays, large peaks and valleys in the sound amplitude 


can be monitored as they progress up the frequency scale. For 
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Tab le V 


Correlation Curves Plotted 


Time Separation Corre lacion 

# of Data Points PeEWeCimeOlnts uw NeneoL Mans Overs | 
50 20 msec 4 Full curve 
" " 4 lst 4% of curves 
ts " 3 Full Curve 
" " 3 IStls Ofpecurves= 
" " 2 Pull eeura.e 
" " 2 2nd 5% of curves 
" L or Full curve* 
" " il lst % of curves* 
95 20 msec 8 Full curve 
" / " 3 lst + of curves 
" . " 3 2nd) 4s On. Curves 
" " 3 lst % of curves 
" " 3 2nd % of curves 
" " 3 3rd % of curves 
" " 3 4th % of curves 


*Data obtained but not plotted. 
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E@eenigher Lan combinations, the incoherency is controlling, 
and in most instances, particular peaks and valleys cannot 
be followed over any appreciable frequency range. 

After an exhaustive examination of the correlation 
Gurves, it was concluded that no information of a significant 
nature could be discerned; it is interesting to note that in 
the 3 fan case a rather rapid de-correlation occurs. The 2 
fan case, in contrast, has a rather high correlation over the 
entire frequency range. This is more or less confirmed by 
visually correlating the curves of the overlays. It should be 
mentioned, that in one instance, the 1 fan combination showed 
slightly less correlation than that of 2 fans. Also of inter- 
est is the fact that the resulting degree of correlation is 
somewhat time-independent as evidenced by the large degree of 
Similarity between the correlations over the full number of 
data points and the correlations over half of the data points. 

In an attempt to quantify the actual effect of fre- 
quency on the scattered sound, graphs of the large, well de- 
fined minima taken from the overlays of Figures 25 to 3l vs. 
time were plotted. The minima values were chosen because, in 
general, they appeared to be more highly defined and localized 
than the rather broad maxima of the curves. The resultant 
minima plots for the various overlays appear in Figures 40 to 
46 and are appropriately identified. Figure 46 illustrates 
how the increasing incoherency makes it impossible to trace 
a logical progression of features as frequency increases. By 


using an assumed extension of the trends developed for the l 
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Figure 40. Minima Plot, 1 Fan, 
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Figure 42. Minima Plot,1 Fan, ? Seconds. 
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observed maxima on the graphical 


x's shown Gn the plot represent 
overlays. 
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and 2 fan cases, and also by assuming that the maxima are, on 
the average, located midway between the plotted minima, it 
was possible to predict the change in frequency required to 
go from a minimum to a maximum at any given instant of time. 
For each of the two curves, this Af was tabulated and aver- 
aged. For the 3 and 4 fan cases, the lack of predictability 
required a more numerical approach. The overlays of Figures 
26 and 29 (3 fans) and the individual graphs for the 4 fan 
cases (not included) were examined and the average frequency 
change from max. to min. determined and tabulated over a num- 
ber of representative cases. The results of both these 
methods are shown graphically in Figure 47. This curve shows 
the predicted change in frequency to go from a minimum to a 


maximum plotted against the average wave height for a given 


fan combination. The 1 and 2 fan frequency changes were 
calculated from 10 values and the 3 and 4 fan cases from 14 

and 17 values respectively. Standard deviations about each 

of these points are also shown in Figure 47. The curve appears 
linear for the smoother surfaces progressing into non-linearity 
by the time the 4 fan surface heights are reached; however 

this lack of linearity of the last point may simply be due to 


a lack of resolution below 10 kHz increments. 
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V. FUTURE RESEARCH 


The author feels that the research described in the fore- 
going pages is a viable means for gaining new knowledge and 
insight into the basic process of scattering of sound off of 
a statistically rough surface. Further investigation should 
provide additional information; a number of approaches to 
improve the previous research are recommended below: 

1. Investigation into the feasibility of using a 5 kHz 
base frequency, instead of the 10 kHz used herein, should be 
conducted. A few trials along these lines were attempted 
with indications that the method is feasible. The frequency 
resolution would then be twice that given in this paper and 
Should be particularly valuable in giving a closer look at 
the 3 and 4 fan surfaces at lower roughness values. Rough- 
ness values less than 1.0 should specifically be looked at. 

a A method for separately looking at the coherent and 
incoherent components of the scattered sound should be de- 
veloped. This appears to be capable of being accomplished 
using the phase information present in the FFT algorithm used 
in the digital computer. The SKIP-1C program has this output 
available. A possibility of merging the capabilities and 
outputs of the SKIP-1C and THCDB(1l) programs should be 
attempted. By removing the incoherent component at the higher 
surface roughness values, more information can be attained 


from the chaos which currently presents itself. 
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3. Alternate means of presenting correlation data should 
be investigated; one possibility would be to calculate the 
correlation of excess pressure amplitudes relative to the mean 
value for a given frequency increment at a given instant in 
time and to display the resultant correlation values as a 
function of time. The advantage of this method is that it 
would lend further support to the hypothesis of consistent 
frequency separation (delta £) between relative maxima and 
minima discussed in section IV.B.5. It is anticipated that 
such curves would reveal a large negative correlation when the 
appropriate frequency of Figure 47 is reached. Correlation 
at twice this frequency increment should result in large posi- 
tive correlation values once again. Analysis of existing 
data arrays in this manner is feasible due to the prepunched 
data cards available in the Ocean Physics Laboratory. 

4, The effects of variation of geometry to include 


larger angles of incidence should be investigated. 
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SKIP-1C: 


MOU : 


oc Llons : 


OueEput ; 


CAL-OlL: 


MOU 


SUEDUtE : 


APPENDIX A: PROGRAM DESCRIPTIONS 


Sample Frequency 

Number of Sample Points per FFT 
Number of FFT's 

Band pass limits desired on printout 
Dee tit lrering 

Triangular hanning function 
Individual power spectral densities 
(PSD) for each FFT 

Normalized average PSD vs. frequency 
over the stipulated bandpass 

Average phase vs. frequency over the 
stipulated bandpass. 

For each FFT, the maximum PSD and 


its associated frequency. 


Number of samples (maximum 8192) 
Input reference voltage set at 

i UGGevoles —( RMS) 

Voltage signal from wavetank probe 


VOUT = AC voltage of wavetank probe. 
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THCDB; 
PApPuUe ; l. Number of sample points per FFT 

2. Total number of FFT's desired (maximum 
Of 50) 

3. Conversion factor for transforming 
PSb aco. VOLES 

4. Amplifier gain 

5. Line loss 

6. Hydrophone sensitivities for every 


Dek? eer OMe” s=lo 0) KHz 


Options: i "SPR or VEroutpute 


Noise run 


SE Put : 1. Mean sound pressure level (or voltage 
level) Vs. frequency (5 kHz increments 
from 0 - 160 kHz) 

2. Standard deviation of SPL vs. frequency 

3. Signal to noise ratio vs. frequency, if 
previous noise run has been performed 

4. Output can be stored on a digital 
cassette or printed on the terminal 


printer 


mHeDB (1): Same as the previous program except the 
program will accept raw A/D input from a 
digital cassette and it can handle up to 


100 FFT's as the associated output. 
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PERK-1A: Program to perform A/D sampling, delay a 
specified time and write the A/D data on 
a digital cassette (to be used as input 


scene VWSEGIRIE (Ike 


Pa put : 1. Total number of A/D samples (maximum of 
8200) 
2. Number of samples per FFT 
3. Sample frequency 
4. Number of FFT's desired 
5. Desired time delay between FFT's 


Output: 1. Raw A/D data stored on digital cassette. 


PRK-MOD: Same as previous program except total 
number of points has been increased to 


a’ maximum value of 13,000 


DORTE ; Progra tOrtakcr Ene sOuLput Of THEDS (1) 
and sort the resulting FFT data by fre- 
quency. slescanwaccepeE sup EO 100 blocks 
of FFT data with up to 75 frequency 


values in each block. 


Output is SPL vs. time for each frequency. 
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